The remarkable warped and twisted gas disk in NGC 3718 
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ABSTRACT 

We have mapped NGC 3718, a nearby bright galaxy in a loose group, and 
its companion NGC 3729 in the 21cm line of neutral hydrogen. NGC 3718 is 
a strikingly unusual galaxy with a strong straight dust lane across the center, 
peculiar diffuse spiral arms, and an extended disk of neutral hydrogen. Earlier 
work showed the gas disk to be strongly twisted, warping through edge-on where 
we see the straight dust lane; stars formed in this gas appear to make up the 'spiral 
arms'. Our improved maps show a twisted but bisymmetric disk of gas extending 
to T or 35 kpc, where the orbital period is roughly 1 Gyr. It is surrounded by 
fragmentary spiral features, and a streamer of gas extending to a cloud lying 12' 
or 60 kpc to the north. We use inspector, a task in gipsy, to fit a tilted-ring 
model interactively to slices through the Hi data cube. The apparent major axis 
swings through 100° from the innermost gas orbits at 30" radius to the outer edge. 
When viewed in the reference frame of the galaxy's stellar disk, the innermost 
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gas orbits are nearly polar, while the outer rings of gas are tilted at 30°-40°. 
The line of nodes, where the gas orbits pass through the plane of the stellar disk, 
twists by roughly 90° about the pole. We do not see gas orbiting in the plane 
of the stellar disk. If we assume that the galaxy's dark halo shares the same 
midplane, then the observed twist can be explained by differential precession in 
a dynamical model in which the dark halo is fairly round. The run of tilt with 
radius is close to what is required for the warped gas disk to precess rigidly in 
the galaxy's gravitational field without changing its shape. This fact probably 
accounts for the longevity of the twisted structure. 

Subject headings: galaxies: individual (NGC 3718) — galaxies: kinematics and 
dynamics 

Facilities: VLA, WIYN, ADS, NED 



Introduction 



The luminous galaxy NGC 3718 (UGC 6524; Arp 214; PRC D-18) and its dwarf 



companion NGC 3729 form a galaxy pair in the loose Ursa Major group (Tully etal. 1996). 



Morphologically NGC 3718 is quite peculiar: see Figure [Tj A strong dark dust lane 



resembling that in Centaurus A (NGC 5128; Dufour etal. 1979) runs almost edge-on and 
straight across the central bulge. Further out, the dust lane diverges into several smooth 
filaments. It then twists by almost 90° into an 'S' shape, forming a diffuse spiral in the 



stellar light which led de Vaucouleurs etal. (1991) to classify the galaxy as a peculiar 



barred spiral. As in Centaurus A, an active nucleus is largely hidden behind the dust lane. 
A compact radio continuum source less than 0.2 pc across has a brightness temperature in 



excess of 10 K (Nagar etal. 2005); Krips etal. (2007) see a jet stretching 0.5" or 40pc to 



the north-west. Optical spectroscopy shows a LINER of Type 1.9 (Ho etal. 1997) with 



weak broad emission at Ha and a strong; narrower line of [Ol] at 6300A. 



The Hi maps of Schwarz ( 1985 ) showed that the dusty gas does not lie in the plane of the 



stellar disk, but forms a complex three-dimensional structure. Some lines of sight pass more 
than once through the gas layer, giving rise to multiple velocity peaks; but the velocity 
field is strikingly bisymmetric. Schwarz was able to describe the gas layer as a violently 
warped disk made up of material following concentric but tilted orbits, which twisted by 
roughly 90° between the inner and outer radii. The strong straight portion of the dust lane 
arises where the orbits turn nearly edge-on to our line of sight, at roughly 200" radius. 



Pott etal. (2004) mapped molecular gas in the inner part of the dust lane using the CO 
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and HCN lines. Krips etal. (2005) combined those data with interferometric maps at ~ 2" 



resolution in the CO 1 — > line at 3 mm. They showed that the CO emission traces an 



inward extension of the warp that Schwarz (1985) derived for the Hi gas. NGC 3718 is 



included as 'related object' D-18 in the Polar Ring Catalog of Whitmore etal. (1990). 



Because polar ring systems contain gas orbiting in more than one plane, these rare objects 
constitute one of the few observational probes of the three-dimensional mass distribution of 



galaxies. Sparke (1990) presented a dynamical model for NGC 3718 to explain the complex 



shape of the twisted Hi disk mapped by Schwarz (1985). According to this model we see 



the underlying disk galaxy almost face-on, with the ring gas in near-polar orbit about it. 
The tilted gas orbits precess about the symmetry axis of the flattened central galaxy and 
its dark halo. Orbits at smaller radius precess more rapidly in the galaxy's gravitational 
field, so the gas disk becomes twisted. This dynamical model reproduced the main features 
of Schwarz's tilted-ring fit. It was consistent with a spherical dark halo, and indicated an 
age for the gas disk of 3 — 4 Gyr. 

To test this model, we mapped the system in the 21 cm line of neutral hydrogen with the 
Very Large Array (VLA) radio telescope. Our new observations, described in Section [2j 



improve on those of Schwarz (1985) in sensitivity and in both velocity and spatial 



resolution. In Section [3] we discuss the gas distribution and compare our results with 
previous observations. In Section [4] we present tilted-ring fits for the gas motions measured 
in Hi and CO, using optical images to resolve ambiguity about where the gas lies in front 
of the stellar body. In Sections [5] and [6] these are interpreted as showing a near-polar disk 
of gas that has become twisted by the differential precession of the gas orbits in the 
galaxy's aspherical gravitational potential. 

Table [T] gives basic information on NGC 3718 and NGC 3729. We adopt the distance of 



17Mpc given by Tully (1998) for the Ursa Major group: there, 1' is equivalent to 4.945 kpc, 

1" = 82.4pc, and 13" = 1.07kpc. NGC 3718 is a very luminous galaxy with 

L B « 3 x 1O 1O L , while for NGC 3729 L B « 7 x 10 9 L Q . |Tully etal] fll996|) classified 



NGC 3718 as T=l (Sa) and NGC 3729 as T=2 (Sab) on the basis of their optical and 
near-infrared images. 



2. Observations in the 21cm Line and Data Reduction 

We used the VLA in the C configuration in four different observing runs in March and 
April 1992, for a total observing time of 26 hours and 28 minutes on source. To obtain the 
required velocity coverage and resolution, we observed using two IFs, tuned at slightly 
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Fig. 1 — Left, B-band image of NGC 3718 taken by E. Wehner with the WIYN 0.9-m 
telescope, showing the strong central dust lane and diffuse 'spiral arms'. North is up and 
East is to the left; the image covers 9' x 13.5'. Right, i?-band image of the central region 
taken by J. Gallagher with the 3.5-m WIYN telescope; the dust lane is close to edge-on, 
with the bright nucleus seen to the north. The dark feature seen closest to the nucleus is at 
PA ss 125°. 
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different frequencies in order to almost double the spectral coverage. Parameters of the 
observations under proposal AS649 are listed in Table |2j Compared to the observations of 
Schwarz (1985), we improve the velocity resolution to 5 kms -1 from 33 kms -1 , and the 



spatial resolution from a 25" x 31" beam to a 13" circular beam. 

The complete data reduction was done using the Astronomical Image Processing System 
(AIPS). The four data sets were calibrated independently, both for amplitude and phase 
gains errors that vary with time, and for those that vary with frequency. The absolute flux 
scale was determined by observing 3C286, which has a well-known flux density. After this, 
the four databases were combined into one. Inspection of a first mapping of the result 
allowed us to determine line-free channels at both edges of the band. The average of these 
channels was subtracted from the uv data set using the AIPS task UVLIN. This new data 
cube, which now contains just line emission, was used throughout in all subsequent 
mapping and cleaning. 

The continuum map shows point sources at the central locations of both NGC 3718 and 
NGC 3729. Table [3] lists positions and fluxes of these point sources, with uncertainties 0.5" 
and 1.0 mJy respectively. The sources coincide to within this accuracy with the positions 



given by Verheijen & Sancisi (2001) at 20 cm and by Krips etal. (2005) at 3 mm; we take 



them to represent the center of each galaxy. The presence of a dust lane prevents an 
accurate optical position for NGC 3718, but the various values in the literature agree with 
this radio position within the error margins. In their Appendix B, |Verheijen &: Sancisi 



(2001) quote 11.4 ± 0.4 mJy for the continuum source in NGC 3718, in reasonable 



agreement with our value of 14.4 mJy; we find 7.9 mJy for NGC 3729, while they give a 
higher flux of 18 ± 0.9 mJy. 

The data were Fourier transformed using the AIPS task IMAGR, using a robustness 



parameter of 0.0, resulting in a circular 13" beam. A proper choice of robustness (Briggs 



1995 Briggs etal. 1999) results in better sensitivity than when using uniform weighting, 



without the strong non-Gaussian beam effects normally associated with natural weighting. 
We made two data cubes: one with the full 5 kms -1 velocity resolution, and once applying 
a smoothing over the frequency axis resulting in 10 kms -1 velocity resolution. The maps 
with 10 kms -1 resolution showed no additional structure in the outer galaxy, so the 
5 kms -1 resolution data were used throughout in this paper. Both cubes were cleaned to a 
1 cr noise level: 0.39 mJy/beam for the full resolution data and 0.30 mJy/beam for the 
frequency-smoothed data, as given in Table |2j 



3. Neutral Hydrogen Datacube for NGC 3718 and NGC 3729 



3.1. Global Results 



We list global Hi properties for NGC 3718 and for NGC 3729 in Table gj For both galaxies, 
we used a method which corrects for the mismatch between the dirty beam and the clean 



beam in the residual map (Jorsater & Van Moorsel 1995) to determine the Hi flux in each 



channel map, and corrected for the attenuation of the VLA primary beam. The global Hi 
profiles in Figure [2] show no sign of absorption against the weak radio continuum source in 
either galaxy. 



For NGC 3718, the Hi flux integral of 118 Jykms 1 agrees with that found by Schwarz 



(1985) and is 20% lower than the value of Verheijen & Sancisi (2001). Our flux integral of 



3.8 Jy kms 1 for NGC 3729 agrees with the latter authors, but is roughly 5 times lower 



than given by Schwarz (1985). A recalculation using the map in Figure 3 of that paper 



yields a much smaller flux integral, so the result quoted in Schwarz ( 1985 ) may have been 
in error. Estimates from single-dish observations vary between 90 Jy kms -1 and 



150 Jy kms 1 (Huchtmeier & Richter 1989). Thus there is little gas in an extended 



component that would be missed in our maps. 

In NGC 3718 we find 8 x 1O 9 M of Hi gas, about twice as much as in the Milky Way, while 
the galaxy is about 50% brighter in stellar light. The ratio Mhi/Lb=0.3, which is about 



average for the sample of gas-rich SO and Sa galaxies studied by Noordermeer etal. (2005). 



NGC 3729 has only 3 x 10 M© of Hi and is gas-poor compared to a normal Sab galaxy; we 



find M H i/L B =0.04 while M ffl /L B =0.1 would be typical (Roberts fe Haynes 1994) 



The regular, steep-sided and symmetric profile of NGC 3718 suggests that the gas has had 
time to settle into a steady state. Between the points at which the emission falls to 20% of 
its peak value we measure a width H / 2 o=476kms~ 1 . If the Hi followed pure circular orbits, 
our measured line width would yield the rotation speed directly: W20 = 2V max sin 2, where 
V max is the maximum rotation speed in the galaxy disk, inclined at angle % to face-on. 
Verheijen & Sancisi (2001) find that in disk galaxies we must subtract about 20 kms -1 



from W20 to correct for random motions in the gas. For NGC 3718 this would imply 
Vmax sin i ~ 230 kms" 1 , with little gas in regions where the circular speed is higher. 

Based on the mean of the velocities at 20% of peak flux, we adopt the systemic velocity 
1/^=995 kms" 1 for NGC 3718 and V sys =W63 kms" 1 for NGC 3729. For NGC 3718, 



Verheijen & Sancisi (2001) derived 993 kms 1 from the midpoint of their Hi global profile, 



and 990 kms 1 by examining the position-velocity diagram along PA=195°. For NGC 3729 
the agreement is even closer: Verheijen & Sancisi (2001) find V,,, s =1060 kms -1 and 
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Fig. 2.— Global profiles of Hi emission in both NGC 3718 (solid line) and NGC 3729 (dashed 
line). Fluxes have been corrected for primary beam attenuation. 
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1063 km s 1 for the two methods respectively. 

3.2. Channel maps 

Figures [3] and [4] show the channel maps for gas in NGC 3718. In a warped disk made up of 
gas on concentric but tilted circular orbits, gas at each velocity above the systemic velocity 
V sys should have a counterpart at the same interval below V sys , at a position point-reflected 
about the galaxy center. In Figure |5j emission from gas in two extreme channels centred at 
765 km s^ 1 and 1222 km s^" 1 has been superposed to show this symmetry. 

The extreme channel maps containing Hi emission are at 755 km s" 1 and 1232 km s -1 , 
separated by almost exactly our measured width W20 = 476 km s -1 . Channel maps at 
5kms _1 lower and higher velocity are empty. The emission in these channels extends from 
30" from the center to 400", so V rot smi should be between 230 km s -1 and 240kms _1 over 
this entire radial range. Within 30", either Hi gas is largely absent, or V ro t sinz is 
considerably lower. 

Within 300" of the center the band of emission in Figure [5] is narrow, suggesting that we 
see the gas orbits within 10° to 20° of edge-on. Tracing the ridge line of the emission in the 
extreme channels then gives us the position angle of the gas orbits, as plotted in Figure [9j 
The kinematic major axis swings from close to PA=100° at 30" from the center to 
PA=190° at radius 400". 



3.3. Blanking and moment maps 

The cube of data can be viewed as a rectangular array of velocity (or frequency) profiles. It 
is standard practice to reduce the spectral line data further by forming maps containing 
the value of the various moments of each profile. The zeroth moment is a map showing the 
spatial distribution of total hydrogen; all velocity information is lost. In calculating the 
zeroth moments we restrict ourselves to that part of the profile where emission is present; 
this avoids contamination of the total HI map with noise. Our method of separating 
emission and noise is automatic: we convolved the cube to a 40" x 40" beam, and masked 
(blanked) all pixels in the high resolution cube which were below a 3<r noise level in the low 
resolution cube. The moment maps are constructed from the unblanked pixels only. This 
method avoids the addition of unrelated noise to the total HI map, and at the same time 
misses little of the low-level HI emission. 

Figure [6] shows the resulting map. Within 100" of the center, we see a narrow dense ridge 
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Fig. 3. — Channel maps for the Hi distribution in NGC 3718 at intervals of roughly 20 km s^ 1 . 
Our adopted systemic velocity V S y S =995kms _1 lies midway between the last channel map 
in this figure and the first in Figure |4| 
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Fig. 4. — Continued from Fig. [3j channel maps for the Hi distribution in NGC 371£ 



- 11 - 




11 33 15 00 32 45 30 15 00 



RIGHT ASCENSION (J2000) 

Fig. 5. — Superposed channel maps for the Hi gas in NGC 3718 at velocities displaced 
roughly 230 km s -1 on either side of the systemic velocity. Lower contours, in red, show gas 
centred at 1221.8 km s -1 ; upper contours, in blue, show gas centred at TGS.Skms^ 1 . 
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RIGHT ASCENSION (J2000) 

Fig. 6. — Total hydrogen in NGC 3718 and NGC 3729, corrected for the primary beam 
attenuation. The lowest contour level is at 3.3 x 10 19 atoms cm -2 or O.26M pc~ 2 , approxi- 
mately at the 3-cr noise level. Higher contours are at 1.0, 3.3, 10.0, and 33.4 10 20 atoms cm" 2 . 
The positions of the central continuum sources are marked with crosses; the beam size is 
indicated by the small circle in the lower left hand corner. 
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of Hi emission along PA=140°, coinciding with the dark dust lane of Figure [T] At larger 
radii this ridge swings counter-clockwise into a 'S' shape. The left panel of Figure [T] shows 
that the diffuse 'spiral arms' that we see in the starlight correspond to regions where the Hi 
density in Figure [6] rises above 3.3 x 10 20 atoms cm -2 or 2.6M pc -2 . This gas density 
barely reaches the threshold of 3-lOM pc -2 normally required for widespread star 



formation in a galaxy disk (e.g. Schaye 2004). The tilted-ring model developed for the Hi 
layer in Section [4] below, and illustrated in Figure [7], implies that the projected density 
along the 'spiral arms' is increased by warping in the gas layer. The true surface density is 
even further below the normal threshold for star formation. 

The stellar 'spiral arms' extend to roughly 250", while the pattern of Hi emission is 
bisymmetric to about 7' or 35kpc from the center, and the Hi disk can be traced to a 
radius of 500" or 41 kpc. To the southeast a spiral-arm fragment extends to a long 
streamer, apparently ending in a gas cloud projected 12' or 59 kpc from the center. This 
arm fragment and a symmetrically placed structure to the northeast are also visible in the 
left panel of Figure [l] as star- forming regions. As in NGC 1058 and NGC 6946 (Ferguson 



etal. 1998; Boomsma etal. 2008 Petric & Rupen 2007), both coherent spiral patterns in 



the gas and continuing star formation are present far beyond the radius where gravitational 
instability should be strong enough to initiate them. 

Smoothing our data cube at 10 kms" 1 velocity resolution further to a 40" beam yields a 
noise level of 0.5mJy/beam. A map of total Hi made from this smoothed cube fails to show 
emission more extended than Figure [6j to a surface density of 0.1 M pc~ 2 or 
1.3 x 10 19 atoms cm -2 . In particular, there is no bridge of emission linking NGC 3718 with 
NGC 3729. 

Beyond 200" from the center, we see only a single velocity peak along each line of sight, 
and the velocity dispersion is generally below 10 kms -1 . Here the first-moment map of 
Figure [7] describes the velocity field of the gas. It shows a pattern characteristic of a 
warped rotating disk: the kinematic major axis (where the velocity is furthest from 
systemic) twists with radius into an 'S '-shape. The orderly rotation and low velocity 
dispersion suggests that the structure is at least a few orbits old. Beyond about 7' or 
35 kpc, the gas of the spiral-arm fragments seen on both sides of the disk along PA=120° 
appears to share in the rotation, although it does not form part of a complete ring. The 
long streamer curving northwards away from the east side of the disk is continuous in both 
position and velocity. Taking the maximum radius as 500" and the circular speed there as 
220 kms -1 (see below) yields a dynamical mass Md yn = 5 x 10 11 M Q , so that M^n/Ljf ~ 7 
in solar units. This is much larger than the value of unity that is typical of an old stellar 



population (Bell etal. 2003); so the galaxy must contain substantial dark matter. 
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Fig. 7. — Mean velocity field (first moment map) of the Hi emission in NGC 3718 and 
NGC 3729; the tilted ring model of Section [4] for the warped gas layer of NGC 3718 is 
superposed. Contours are spaced at intervals of SOkms^ 1 around the systemic velocity 
of 995 km s -1 . All the contours in the gas streamer on the northeast side of the disk of 
NGC 3718 are at 795 km s -1 . Beyond about 200" from the center, where velocity profiles are 
singly peaked, the values in this map are representative of the radial velocity of the gas at 
that position. The beam size is indicated by the small circle in the lower left hand corner. 



- 15 - 



The small galaxy NGC 3729 is projected 11' to the east, and shows a clear signature of 
rotation in gas that extends to 1' or 5kpc radius. ^From their K-band images, |Tully etal. 



(1996) find an isophotal ellipticity e = 1 — b/a = 0.32. If we assume the disk to be round 
(although the galaxy is classified as barred), it is inclined 48° from face-on. Further 
assuming that the Hi gas shares this plane yields a dynamical mass M-dyn = 35 x 10 9 M Q 
and M/Lk ~ 2 in solar units. Here we cannot draw strong conclusions about the presence 
of dark matter. 

In the central parts of NGC 3718 the velocity profiles are neither singly-peaked nor 
symmetric. Figure [8] shows a position- velocity plot through the center at a position angle 
of 140°, along the ridge of bright HI emission visible in Figure [6] This plot is highly 
symmetrical, as expected for gas in circular orbit about the galaxy center. There are two 
main components: the very strong inner one shows velocities rising steeply to 230 km s -1 at 
80" from the center, while in the outer component rotation speeds increase almost linearly 
to 60kms _1 at 300" radius. We interpret the slower-rotating gas as following an orbit at 
larger radius; we see only the portion projected close to the center, where the radial 
velocity is small. On the western side and at negative velocities, there is a third and much 
weaker component with an intermediate slope. Looking towards any point along this line 
within 80" of the center, we would see a double or triple peak in the Hi velocity profile. If 
the Hi gas forms a continuously warped disk and we look once through it in the outer 
parts, then each line of sight must cross the gas layer an odd number of times, so we expect 
triple profiles. There is a slight indication that the weak third component may have a 
counterpart to the east at positive velocities. 



Krips etal. (2005) made interferometric maps of the molecular gas associated with the 



inner part of the dust lane. They combined several pointings with the single-dish 



observations of Pott etal. (2004) to probe the dust lane to 70" radius with ~ 2" resolution 
in the CO 1 — >• line at 3 mm. Their Figure 17 displays a position-velocity diagram along 
an axis at PA=130°. Like our Figure [8j it is highly symmetric about the center, with 
emission at velocities rising to 220 km s^ 1 at 70" radius, and an even smaller-scale structure 
where speeds reach 250 km s -1 within 10" of the center. This nuclear component would 
correspond to an edge-on disk of diameter 1.5 kpc. If atomic and molecular gas share the 
same kinematics, then Hi must be largely absent within 30" of the center; otherwise, our 



Figure [5] would not show a gap in high- velocity emission close to the center. Reshetnikov 



Combes| ( 1994 ) measured velocities in the Ha line of ionized gas in the central regions of 
NGC 3718. Along PA=130°, they found velocities rising to 260 ± 20kms" 1 within 80" of 
the center, which is consistent with the results in CO. 
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Fig. 8. — Position-velocity plot through the center of NGC 3718, and along the apparent HI 
ridge at PA=140°. The axes are labeled relative to the systemic velocity V sys =995kms~ 1 
and the center of NGC 3718: right ascension increases to the left. The lowest contour is at 
0.75 mJy/beam; higher contours are odd multiples. The sense of the velocity axis is chosen 
for comparison with the figures of Krips etal. (2005). 



4. Tilted-ring models for the Hi gas 



4.1. Fitting a tilted-ring model 



Because of the very high degree of symmetry in the channel maps of Figures [3] and [4], we 



follow Schwarz (1985) in modeling the gas within a radius of 500" as a strongly- warped but 



otherwise symmetric disk. The disk is made up of rings of material, following near-circular 
orbits that are concentric but tilted. Because the emission does not peak symmetrically 
about a mean velocity at each point on the sky, we cannot use tasks such as ROTCUR 



(Begeman 1987, 1989) to determine the ring orientations by fitting to the mean velocity 



field, as measured by the first-moment map. Instead, we built the task inspector in 



GIPSY (Vogelaar & Terlouw 2001) to compare the predictions of such a model to various 



two-dimensional cuts through the three-dimensional cube of data. 



Following the convention of Rots (1975) and Begeman (1989), we measure the position 



angle p of each gas orbit anti-clockwise from north to the line of nodes (the kinematic 
major axis) on the receding side of the galaxy. Note that this definition of p is 180° 



different from that of Schwarz (1985). The orbital inclination i runs from zero as the spin 



axis points towards the observer, through 90° for an edge-on ring, to 180°. 

Neglecting the effect of both random motion in the gas and our finite beam size, 
inspector calculates the expected velocities at which each ring of Hi should contribute to 
a given longitude- velocity cut, or the positions at which its emission should appear in a 
given channel map. Fixing the central velocity at 990 km s -1 gave a slightly better fit than 
the central value of 995kms _1 that we derived from the global profile. We placed the ring 
centers at the radio continuum source, and adjusted the rotation speeds and the ring angles 
interactively, using inspector to compare model predictions with the position-velocity 
cuts and channel maps. 

Figure [9] shows that the position angle of the gas orbits is the best-determined quantity. In 
the top panel, we see that at radii r > 40" the results from inspector are in excellent 
agreement with those obtained by tracing the ridge line of emission in the extreme velocity 



channels. They match fairly well to the model fit by Krips etal. (2005) to the CO 



observations: see below. Near the galaxy center, the apparent major axis of the gas orbits 
lies roughly east-west, almost orthogonal to the major axis of the stellar light. It then 
turns counterclockwise towards north-south, twisting quite sharply at smaller radii and 
then more slowly beyond 300". 



The most difficult quantity to determine is the rotation speed. Initially, we used the 
rotation curve fit by Schwarz ( 1985 ) to the earlier Hi observations. This yielded the model 
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Fig. 9. — Tilted ring model for the warped Hi disk in NGC 3718. In the top panel, red 
lines and open circles show the position angle p of the receding line of nodes, as defined 
by the ridge lines of intensity along the three extreme channel maps above and below the 
systemic velocity (at 754.9 km s" 1 , 760.1 km s" 1 , 770.5 kms -1 , 1232.2kms" 1 , 1227.0 kins" 1 
and 1216.6 kms -1 ). Filled circles and a dashed or solid line show our tilted-ring models 
inspector! and inspector.2 respectively. The crosses and green short-dashed line shows 



the tilted-ring fit derived by Schwarz (1985), while the triangles and the blue (dotted) line 



show the fit from Figure 6 of Jozsa etal. (2004). The purple long-dashed line shows the 



warped disk of CO from Krips etal. (2005). The second panel shows inclination i, and the 



third panel the assumed or fitted circular speed V rot . The bottom panel gives V rot smi, the 
maximum speed along the line of sight. The red lines with long and short dashes show the 



constraint derived in Section 3.2 above above, and the red star shows the innermost velocities 



seen in CO by Krips etal. (2005). 
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inspectorI. However, the CO velocities of Krips etal. (2005) show a rise to 250 km s 



within 10" of the center. Our 13" synthesized beam is roughly 1 kpc across, so we do not 
expect to resolve a rapid central rise in the rotation curve. For the model INSPECTOR^, we 
began our iteration with V ro t set at 250kms _1 , and decreased it in the outer parts only 
when we could not otherwise obtain a good fit. Figure [7] shows the geometry of this 



tilted-ring model; the run of position angle and inclination are given in . Figures 10 - 12 
compare the model predictions with position-velocity cuts and channel maps. 

The lower panels of Figure [9] show the runs of inclination and rotation speed. The 
multiply-peaked velocity profiles illustrated in Figure [8] require that the warped gas disk 
passes through edge-on with i = 90°; we place this transition between 160" and 220". 
Within this radius, lines of sight can pass three times through the disk. The position angle 
in this region of nearly edge-on gas orbits is 155°-175°, running along the dark central dust 
lane in the left panel of Figure [TJ The velocity field of the Hi gas is exactly the same for a 
ring of inclination i and one at 180° — i; we use the dust distribution in the right panel of 
Figure [T] to resolve this ambiguity. There, we see the bright nucleus to the north of the 
dust lane; so the south side of the dusty gas disk is closest to us. The gas recedes on the 
east side, so its spin axis points away from us, meaning that i > 90°. The warp appears 



smooth as the disk twists through edge-on, so we follow Schwarz (1985) in assuming that 



the inclination decreases monotonically to i < 90° at larger radii. 



Schwarz ( 1985 ) constrained the position angle of the edge-on gas orbit by examining the 
emission peak closest to the center at velocities close to systemic. He found that the 
centroid of that peak moved along a line in PA = —23° as the velocity decreased through 
V sys . This is the behavior expected along an edge-on circular orbit in PA = 157°. We 
repeated this exercise for the present data set as a consistency check, and find that the 
central peak moves along PA = 135°. Our measured velocity gradient corresponds to a 
ring at radius 90", where the top panel of Figure [9] shows that the position angle indeed 
reaches 135°. Thus the gas orbit at this radius is very close to edge-on. 

Within 300" of the center the Hi gas orbits are less than 10° from edge-on. So the 
measured speed V rot sin i should be very near to the orbital speed itself. Our stipulation 
that V ro t should decrease monotonically means that our model predicts too high a value for 
V ro t sin i. To avoid this we would need a rotation curve like that of the model inspector!, 
which peaks beyond the optical radius R25. Closed loops in the channel maps at 
1190.6 km s -1 in Figure [3] and at 796.4 km s _1 in Figure [4] show that either the rotation 
speed must drop in the outer disk, or the gas orbits turn closer to face-on. We find that 
both effects are present. At large radii the shape of the total HI map in Figure [6] shows 
that the outermost gas orbits turn to i ~ 60°, or ~ 30° from edge-on. They cannot become 



-20- 



much more face-on, or the predicted east-west extent of gas in channel maps near the 
systemic velocity becomes much larger than observed, and the two arms of the "fork" in 



Figures [10] and [TT] are too wide-open. To reproduce the closed loops, we had to reduce the 
model rotation speed to about 220 km s -1 near the outer edge. This behavior is consistent 



with the 10%-30% drop in rotation speed that Noordermeer etal. (2007) found to be 
common in massive SO and Sa galaxies, with rotation speeds above 200 km s -1 . 

Far from the galaxy center, emission in the channel map near 990 km s -1 extends almost 
east-west for about 200" on each side of the center. However, the fit at this velocity (shown 



in the top left panel in both Figures 10 and 11) is better on the east (left) side than the 



west. Also, Figure [6] shows that the gas furthest to the east and west does not seem to be 
part of a complete ring. So we treat our model with caution within 40" and beyond 400" 
radius. 

We see in Figure [9] that the molecular gas follows the same warped disk structure as the 



inner Hi layer. Krips etal. (2005) fitted a tilted-ring model to describe the CO kinematics. 



They chose a model rotation curve close to that of Schwarz (1985): the rotation speed V(r) 
is taken as 235kms _1 within 100" of the center, rising linearly to 255 km s" 1 at 130". Their 
Figure 12 displays the derived run of tilt angle with radius, relative to a reference plane 
inclined by 70° (or 110°) to the plane of the sky, and with the approaching line of nodes at 
PA=-60°. With respect to that reference plane, their model takes the twist angle to 
increase with radius r as twist oc cos(tilt) x r/V(r) (compare Equation [I] below). Taking 
the reference inclination as % = 110° and using the tilt and twist angles kindly supplied by 
Dr. Krips, we recovered the inclination and position angles of their model relative to the 
sky plane, as shown in Figure [9j The position angle agrees well with what we derive from 
the Hi observations. The inclination oscillates because of the form that they chose for the 
twist, but the product V rot sin 2 is very close to that for the Hi layer. 

These sets of derived quantities each represent an eyeball fit to a constrained parametric 
model. This contrasts with systematic fitting techniques such as ROTCUR that are applied 
to galaxies with a single-valued velocity field. Because the velocity field represents an 
integral over the full data cube, it is smooth and relatively insensitive to the clumpy 
distribution of emitting gas, and can be compared directly with a model in which gas 
orbits are uniformly filled. In a galaxy like NGC 3718 we must work with the full 3-D data 
cube, where the patchy emitting gas lies close to a warped and folded 2-D surface. 



Jozsa etal. (2004) present a model for the Hi layer in NGC 3718 from TiRiFiC, a new 



method (Jozsa etal. 2007) which fits a tilted-ring model automatically to the full data 



cube. Their observations at Westerbork had a resolution of 12". Figure [9] shows that the 
run of position angle is very similar to ours, and the inclination shows the same decreasing 
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Fig. 10. — Tilted ring model inspector^ compared with channel maps at velocities greater 
than our adopted central velocity of 990 km s -1 . Crosses indicate emission from each of the 
model rings; the size of the cross increases proportionally with the ring radius. The central 
velocity for each map is given in kms" 1 in the top left corner; maps with the larger labels are 
displaced from the central velocity by the same amount as the maps in corresponding panels 
of Figure [IT] Other channels are chosen to illustrate features such as the closed contours at 
l^kHis" 1 . 
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Fig. 11. — As Figure 10, but for channel maps at velocities below the central velocity. Maps 
with the larger velocity labels are displaced from 990 km s" 1 by the same amount as the 



maps in corresponding panels of Figure 10 
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Fig. 12. — Tilted ring model compared with longitude-velocity cuts taken along the east-west 
direction. The velocity of the front (closest) portion of each ring is shown by an open blue 
square, and the rear (more distant) portion by a red cross. Large squares enclose symbols 
at radius 400"; large triangles show rings at 300". Large red and blue circles with inscribed 
crosses show rings at 200"; in the central cut only, these circles almost coincide and are 
shown as a single light symbol, very close to the central velocity. Rings at 100" appear in 
the central cut only, as large circles with central dots. 
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trend within 400". The run of inclination differs; Jozsa etal. (2004) estimate that the gas 
orbits turn through edge-on closer to the center, at 80" radius and in PA w 130°, and that 
at larger radii the gas remains further from edge-on than indicated by our model-fits. In 
the central 240" the run of V ro t sin i falls below the constraint that we derived from the 
channel maps of Section 3.2. The implied rotation curve is not monotonic, rising from 
210 km s" 1 with maxima at 120" and 320". The differences between the sets of curves in 
Figure [9] illustrate the difficulties of the fitting methods, the limitations of the data, and 
deviations from uniformly filled concentric circular orbits. 



4.2. Relation between the gas layer and the stellar disk 



In Figure [TJ we seem to see the stellar disk of NGC 3718 close to face-on. We cannot easily 
measure its orientation from the isophotes, because of the obscuring dust. From 



near-infrared photometry in the H band (1.6 /xm) within 50" of the center, Peletier & 
Willner (1993) find isophotes elongated in PA=112°, with ellipticity e = 1 — b/a — 0.17 (see 
their Table 4). |Tully etal] (|1996[) give e = 0.58 at PA=195°, measured between 150" and 



250" radius (see their Figure 8 and Table 2), which would correspond to a round disk seen 
55° to face-on (assuming an intrinsic axis ratio b/a = 0.2). In fact the elongation is caused 
by the stellar light of the 'spiral arms' seen in Figure [T] Marc Verheijen kindly supplied us 



with two images at K band (2.2 /xm) taken by Tully etal. (1996); only the lower-reslution 
image with 2.052" pixels was used in their paper. Neither shows any sign of the dust lane, 
even at the center. ^From their higher- resolution image with 0.753" pixels, we measure an 
ellipticity e = 1 — b/a = 0.11 to 0.12 at 25"- 27", which is within the first exponential scale 
length of the disk but beyond most of the bulge light (see below). This would correspond 
to a round disk with intrinsic b/a = 0.2 inclined 28° to face-on. The major axis at 
PA « 12° is almost the same as at large radii. However, oval distortions of 10% are 



common in galaxy disks (Rix & Zaritsky 1995 Kormendy & Kennicutt 2004), especially 



among earlier types (Ryden 2006). In what follows, we assume that the position angle 
where the galaxy disk intersects the plane of the sky lies in the range p g = 195° ± 20 °. 



Heraudeau & Simien (1998) measured velocities along PA=15°, and find a rise to roughly 



100 km s -1 at 20"-30" radius on both sides of the center. According to their Figure 1 the 
southwest side of the stellar disk is receding, just as for the outer Hi, so the receding line of 
nodes lies near PA=195°. If the stellar disk is indeed inclined at i=28°, then for circular 
speeds close to 250 km s" 1 we would expect to see motions of about 110 km s -1 along the 
kinematic major axis. So these observations are consistent with a round stellar disk 
inclined with its apparent major axis close to their slit position. |Heraudeau fc Simien 
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(1998) find a central velocity dispersion of 193 km s 1 (including the factor f bulge i n their 
Table 1), which drops to about 100 kms -1 at 30". 

From this photometric and kinematic evidence, the galaxy disk appears to be nearly 



face-on, as suggested by the dynamical model developed by Sparke ( 1990 ) for the warped 



and twisted gas layer. However, that model took the plane of the stellar disk to be inclined 
with i n 



Sparke 2002). Over most of its 



u g ~ 20° at a position angle close to PA=-90° (see 
radial extent, the Hi disk is then tilted by about 80° with respect to this reference plane, 
and its twisting could be explained by differential precession. But if the stellar disk indeed 
has this orientation, we would expect streaming speeds to be low along the direction 



PA=15° explored by Heraudeau & Simien (1998). It seems more likely that the stellar disk 



intersects the sky plane along a line closer to PA=15° (or equivalently PA=195°). 
Accordingly, we abandon the earlier model. 

The gas orbits of our tilted-ring model nowhere lie close to the plane of the galaxy's stellar 
disk. The disk of NGC 3718 seems to be that of an SO galaxy, substantially free of cool gas. 



This is similar to NGC 2655 (Sparke etal. 2008), an SO/a galaxy with a strong asymmetric 
central dust lane. It contrasts with NGC 660 (PRC C-13), a starburst galaxy with a 



twisted polar ring that is tilted by roughly 55° to the stellar disk (van Driel etal. 1995). In 



NGC 660 the host galaxy's disk is as gas-rich as a typical Sc galaxy, and contains a quarter 
of the Hi gas in the system. 



Although it contains 8 x 10 M of Hi gas, with 4 x 10 M© of molecular material (Pott 



etal. 2004), the modest far-infrared luminosity of 5 x 10 L Q (Rice etal. 1988, with our 



adopted distance) shows that this galaxy is making few new stars. The compact star 
clusters visible in the right panel of Figure [l] are bluer than their surroundings ( Trinh et al. 



2006), indicating relative youth, and the stellar 'spiral arms' noted in Section 3.3 show that 



some starbirth still occurs in this galaxy. But because the cool gas is not concentrated 
towards the main stellar disk, its density may be too low for efficient star formation. 



5. Why should the gas layer be warped and twisted? 

What might have caused the gas layer in NGC 3718 to become warped and twisted? A 
disk of material following orbits tilted away from the galaxy equator will tend to twist 
because of differential precession. In an oblate galaxy, consider a cloud of gas following an 
orbit tipped by an angle a away from the equator, passing upward through the midplane. 
The cloud will make a complete vertical oscillation and again cross the midplane traveling 
upward, before it has made a whole orbit about the center. The tilt of its orbit remains 
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constant, but the line of nodes, where that orbit crosses the symmetry plane, regresses in 
the direction opposite to the orbital motion: (see e.g. Section 5.8 of Goldstein etal. 2002). 



The angular precession rate Q p for an orbit at radius r inclined by an angle a is related to 
the circular speed V(r) by 

1 



a 



rV{r) d cos a 



-e$ cosay(r)/r or — e$cosaf2(r) 



1) 



Here ($) represents the gravitational potential energy, averaged over the ring {e.g. Sparke 
1986), and fl(r) = V/r is the orbital angular speed. The quantity e$ measures the 



flattening of the potential; it is positive for an oblate system, so Q p is negative. Because 
the orbital periods are shorter towards the center the inner orbits will regress faster, unless 
the galaxy's flattening increases strongly with radius. Thus a gas disk made up of material 
on concentric tilted orbits generally develops a leading twist. Conversely, a disk in a 
prolate galaxy potential will twist in a sense that trails the rotation. 

Its very regular velocity field suggests that the outer Hi disk of NGC 3718 has been in 
place for at least a couple of orbits. For a rotation speed of 230 kms" 1 , the orbital period 
at 400" is roughly 900 Myr, implying an age of at least 2 Gyr. Rotation times in the inner 
disk are much shorter, and at 40" radius this would correspond to at least 20 orbits. The 
position angle of the gas orbits has twisted by about 120° between these radii. If that twist 
represents precession in a system of roughly constant flattening e$, then by Equation [T] we 
must have \Q P \ < Q/60 for the inner orbits, or e$cos« < 1/60. At first glance, this implies 
that the gravitational potential must be improbably spherical to prevent the disk from 
twisting around itself many times in its > 2 Gyr lifetime. 

The twisted dust lane of the peculiar SO galaxy NGC 4753 presents a similar dynamical 



problem. Here Steiman-Cameron et al. ( 1992 ) found a gas disk extending to roughly 7 



times the radius of the inner edge at 13" (lkpc at an assumed distance of 15.8 Mpc) that 
appears to wrap by almost two complete turns around the galaxy pole. They argue that 
the outer disk is at least six orbits old, corresponding to 40 orbits at the inner edge. 
Precessional twisting then implies a nearly spherical mass distribution with axis ratio 
b/a> 0.84. The stellar body is flattened with an axis ratio roughly 2:1 (b/a = 0.5), so 
these authors conclude that the dark halo must be both round and gravitationally 
dominant even well within the main stellar body of this luminous galaxy. 



A different model for the warped disk of dusty gas in Centaurus A was proposed by |van 



Albada etal. (1982): this had the great advantage of representing a stable equilibrium 



state. It requires that the galaxy's mass distribution is not axisymmetric, but a triaxial 
spheroid tumbling about its short axis. The potential then supports a family of stable 



anomalous orbits' described by Heisler etal. (1982), which make up a warped disk. The 
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dusty gas should settle onto these orbits, in a process studied numerically by |Habe"S 



Ikeuchi (1985), Habe & Ikeuchi (1988), Steiman-Cameron h Durisen (1988) and Colley & 



Sparke (1996). Near the galaxy center, in the core of the gravitational potential, the 



anomalous orbits circle the long axis. The pole tilts with radius, until at large radii the 
orbits lie in the 'equatorial' plane perpendicular to the short axis, circling it in the opposite 
sense to the one in which the figure tumbles. The anomalous orbits make up a twisted disk 
which follows a restricted warp: orbits at all radii cross a single line of nodes, which at each 
instant lies along the intermediate axis of the tumbling triaxial galaxy. 

We can understand the anomalous family as a set of orbits that precess about the short 
axis of the galaxy at exactly the right rate to keep up with the tumbing galaxy potential. 
Orbits in a triaxial galaxy will precess about the long or short axis at an average rate given 
by Equation [TJ where e$ is now the average flattening about that axis {e.g. 



Steiman-Cameron & Durisen 1984). When the triaxial figure tumbles about its own short 



axis at a rate f^, the stable anomalous orbit family consists of just those orbits that 
precess at the rate Q p = Q t . For orbits circling the short axis we have e$ > 0, so 
Equation [T] requires the orbital motion to be retrograde with Q t < 0. If the system is 
equally aspherical at all radii and V(r) is constant, then Q p is constant when 

cos a oc r . 



(2) 



As Heisler etal. (1982) point out, the anomalous family tilts over to reach the galaxy's 



equatorial plane at the radius where the rate Vt p of free precession for an orbit that is only 
slightly tilted from the equator becomes equal to the tumbling speed fl t - 

To test whether the anomalous orbit family can represent the warped Hi layer of 
NGC 3718, we specify the position angle of each gas orbit by the unit vector 1 along the 
receding line of nodes, where the orbit crosses the plane of the sky. The spin axis is along n 
which is perpendicular to 1, and we take m in the plane of the ring to complete the 
right-handed set 1, m, n. We take Cartesian coordinates x, y, z from the galaxy center, 
with z pointing towards the observer, x to the east and y to the north. In these coordinates, 
the vectors 1, m, n are related to the inclination i and position angle p of Section [4] by 

1 = (— ship, cosp, 0) , m = (— cos 2 cosp, — cosi ship, — sini) 



and n = (— sin i cosp, — sin i ship, cos i) 



(3) 



The planes defined by two rings with normals along vectors ri! , n 2 intersect along the line 
n i x n 2- A restricted warp is one in which this vector points in the same direction for all 
pairs of rings. 

We can specify the equatorial plane of the galaxy by the apparent position angle p g and 
inclination i g of a circular ring lying in that plane. (Note that if the stellar body is triaxial, 
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the position angle of the galaxy's apparent major axis may differ from p g .) Defining the 
corresponding vectors \ g , m g and n g , the angle tilt between a gas ring and that plane is 
given by 

cos(tilt) = n • n g = cos i cos i g + sin i sin i g cos(p — p g ) . (4) 

The ring intersects the galaxy's equatorial plane along the direction n x n g . We define the 
twist to be the angle in the equatorial plane between n x n g and the vector \ g where the 
equator intersects the sky plane; so 



sm(tilt) cos(twist) = l g ■ n x n g = [cosi g sin 2 cos(p — p g ) — cosi smi g ] = n ■ m 



■9 ■ 



and sm(tilt) sm(twist) = m g ■ n x n g = sini sin(p — p g ) = — n • 1 



9 ■ 



(5) 



With these definitions, tilt = i and twist = p — p g when the galaxy's equatorial plane 
coincides with the plane of the sky so that i g = 0. These are related to the angles 6,(3 oi 
Schwarz (1985) by tilt = 6 and twist = 180° — (3. The pair of angles (tilt, twist) describes 
the same ring as (—tilt, twist + 180°). Just as for the inclination i, we usually take 
< tilt < 180° so that sin(tilt) is positive. 



The solid curves of Figure 13 show the angles for the Hi orbits, relative to a plane with 



=95° and p g =105°, close to that of the gas ring at r = 40", our innermost 



reliably-determined orbit. As Schwarz (1985) found, over the best-measured portion of the 



disk the rings fall close to a restricted warp. (Schwarz's reference plane corresponds to 
i g = 104°, p g = 114° which is close to the position angle of our gas at 50" radius.) The twist 
angles of all the rings between 40" and 300" fall within 10° of a common value. If we take 
the ring at 40" to define the polar plane of the galaxy's potential, then in 40" < r < 400" 
the orientation of the gas orbits changes almost exactly from polar to equatorial, as 



predicted by the model of van Albada etal. (1982). 



However, there are two difficulties with interpreting the gas motions as material following 
anomalous retrograde orbits. First, the outer Hi orbits should lie perpendicular to the 
short axis of the triaxial potential. Simulations combining dissipative gas with cold dark 



matter (Dubinski 1994 Kazantzidis etal. 2004 Bailin etal. 2005) find that this axis tends 



to be perpendicular to the galaxy disk. Those models are directly applicable to luminous 
galaxies like NGC 3718, where the stellar disk should dominate the gravitational force 



within the optical radius (e.g. Kassin etal. 2006). But we saw in Section [4] that the 
luminous disk is close to face-on, unlike the outer Hi orbits. Also, the shape of the warp 
does not follow the prediction of Equation [2j given by the straight dashed line in the 
bottom panel of Figure 13 As Schwarz ( 1985[ ) noted, within 200" of the center the tilt of 
the disk changes too rapidly to fit this description. 



The stellar disk appears to be tipped by about 28° from face-on (Section [i] above). So we 
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Fig. 13. — For the ring model inspector^, angles with respect to a reference plane at 
inclination i g and position angle p g . The top panel shows twist measured relative to the ring 
at 140". The red solid line refers to the 'restricted warp' obtained for i g =95° and p g =105°: 
the twist is nearly constant in the range 100" < r < 400". The blue line with dots refers 
to the orientation of the stellar disk implied by the K-band isophotes: i g =28°, p g =195°; 
the twist is leading. The green dash-dotted line is for i g =152° and p g =195°, the other 
possible orientation corresponding to the K-band isophotes. The twist now has a trailing 
sense. The middle panel shows the angle a between each ring and the galaxy equator in 
the corresponding model: a = 90° — tilt for the restricted warp and a = tilt for the other 
models. The ring inclination decreases monotonically with radius from polar to equatorial 
for all the models except the last. The bottom panel shows cos a; the dashed line shows the 
relation cos a oc r of Equation [2] The run of twist and tilt for the 'restricted warp' model 
and that with i g =28°, p g =195° is given in Table [5] 
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have either i„=28° or i g —152°, depending on whether the east or the west side of the disk 



is closer to us. The blue curves with dots in Figure 13 shows that for the combination 



i g =28° } p g =195° } the central gas disk is very nearly polar while orbits at larger radius tilt 
monotonically towards the galaxy plane. The bottom panel shows that it follows rather 
closely the curve cos a oc r that we expect for the anomalous orbit family. However, this is 
far from a restricted warp: the twist angle increases by about 120° between the inner and 
outer radii. The twist has a leading sense relative to the orbital motion, as we expect for 
differential precession in an oblate galaxy potential. 



When we choose z„=152°, the dash-dot line in Figure 13 shows that the tilt is not 



monotonic. The gas orbits are nearly polar near the center, then dip by about 20°, and 
then warp up towards the pole and over it at 420". Because the flattened stellar body of 
the galaxy should dominates the gravitational force within 150" of the center (see below), 
the potential should be oblate and we expect the precessional twist to have a leading sense. 
Instead, we see a trailing twist. The shape of the gas layer is neither a stable configuration, 
nor a natural result of precessional twisting. We do not consider this model further, but 
adopt i g =28° for the stellar disk. 



If the ring is twisted about the galaxy pole, it cannot be in a steady state: e.g. |Hunter & 



Toomre (1969), Sparke (1986), Arnaboldi & Sparke (1994). Instead, the gas orbits suffer 



differential precession according to Equation [T] In the following section, we construct a 
mass model for the galaxy, to examine how fast the gas layer should twist up, and for how 
long the warped gas disk might have been in place. This model is similar to those presented 



for Centaurus A by Quillen etal. (1992), Quillen etal. (1993) and Sparke (1996), where the 



complex warped structure results from an interplay between self-gravity and precession. 



6. An illustrative dynamical model 



We now examine how a tilted gas disk would precess in a simple axisymmetric mass model 
for the disk, bulge and dark halo of NGC 3718. Our model for the stellar component is 
based on near-infrared photometry, to minimize the effect of dust absorption. From deep 



K-band images that trace the galaxy's light beyond 300" from the center, Tully etal. 
(1996) measure a scale length Hr = 56.6" = 4.66 kpc. This is longer than the Hr = 27" 



found by Peletier & Willner (1993) in the H band, and by Chitre & Jog (2002) from 



2MASS K-band photometry, but both of these images were much shallower. Making our 



own ellipse fits to the published image of Tully etal. (1996) with 2.052" pixels confirms the 
longer scale length; so for our illustrative model we adopt h R = 55". We follow fTully etal. 



(1996) in taking the stellar disk to have an intrinsic axis ratio b/a = 0.2, and calculate the 
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forces from this thickened exponential disk as described in Sackett & Sparke ( 1990 ) 



We take the bulge to be spherical. The K-band radial profile in Figure 8 of |Tully etal. 



(1996) appears roughly exponential outside 20", as does the 2MASS profile measured by 



Chitre & Jog (2002). Since our innermost measured Hi orbits lie further out, it does not 



matter how we distribute the bulge mass within that radius. For simplicity we model the 
bulee Plummer sphere with core radius rp = 10". 

The rotation curve of Figure [9] remains nearly flat to 400", which is at least four scale 
lengths of the stellar disk. This, and the high mass-to-light ratio M.dyn/'^K ~ 7 
(Section [33), requires an extended dark halo. We use the pseudo-isothermal form of 
Sackett & Sparke (1990), parametrized by the flattening e = 1 — b/a of the equidensity 



contours, the core radius r# and the asympototic circular speed Vh- For a given halo 
flattening, we set r# and Vh by requiring that the combined rotation curve V(r) from the 
bulge, disk and halo remains approximately flat. To calculate V(r) we use the equatorial 



rotation curve of the halo from Equation 4 of Sackett & Sparke (1990), but we average the 



inward pull of the exponential disk over a circular ring at the appropriate tilt angle. The 



halo torque is computed as described in Sackett & Sparke (1990). The torques from the 



halo and the flattened disk are added to calculate the precession rate according to 
Equation [1} We do not include the HI gas mass in calculating the rotation curve: see 



below. Models like that of Sparke (1996) for the warped disk in Centaurus A shows that 



the self-gravity of the warped disk can also affect details of how it resists precessional 
twisting. In this case the disk is very strongly twisted, so this effect is likely to be small, 
and we do not include it. 



The top panel of Figure 14 shows the rotation curve from this model. The disk mass 
= 5 x 1O 1O M and the halo is chosen to have a small core radius, r#=10", so the 
rotation speed declines gently with radius. To provide a flat rotation curve at the center, 
the bulge mass is relatively small, M& = 2 x 10 10 M Q ; the overall mass-to-light ratio 
M/Lk = 1 in solar units. This is a 'maximim disk' model: the disk and bulge must 
dominate the rotation curve within 2Hr to provide the observed declining rotation curve. 
We initially take the dark halo to be spherical. 

The middle and bottom panels show the orientation of the gas orbits in our fit 



inspector^, relative to a 'galaxy' oriented with i g = 28° and 



195°. As discussed in 



Section [4] above, the viewing angles i g and p g for the stellar disk are also uncertain. The 
predicted twist is not very sensitive to a change in the inclination i g , but decreasing the 
position angle p g slightly will change the sign of cos(tilt) for the inner, near-polar orbits, 
and hence the sense of their precession. The open circles show the orientation when we 
take p„=175°. 
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Fig. 14. — Rotation curve and expected twisting for a dynamical model with = 5x 1O 1O M 
and Mb = 2 x 1O 1O M . The spherical halo has r#=10" and Vff=200kms _1 . Top: points 
show the rotation curve of our tilted-ring fit INSPECTORS in Figure [9j curves show the 
total rotation predicted from the dynamical model (solid), with the contributions of dark 
halo (dotted), disk (dashed) and bulge (dash-dot). Middle: angle tilt of the Hi orbits from 
inspectors relative to a stellar disk with i g =28°, p 9 =195° (dashed line with filled dots) and 
for p g =175° (open circles). An exactly polar orbit has tilt=90°. Bottom: angle twist for the 
Hi orbits (filled dots for p g = 195°, open circles for p g =175°), and the precessional twisting 
predicted by the mass model (line with dots for p g =195°, line with crosses for p g =175°). The 
dashed curve shows the result for p g = 195°, when the halo is flattened with an E3 shape. All 
twists are measured relative to the gas orbit at r = 240". Measured angles and consequently 
the predicted twists are uncertain within 40" radius. 
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The bottom panel shows how much twisting a gas disk would suffer over our inferred 
minimum 2 Gyr lifetime, if it was initially warped but not twisted, so that its tilted orbits 
intersected the 'galaxy' plane along a single straight line of nodes. Comparing this 
prediction to the twist angles derived from our tilted-ring fit, we see that after 2 Gyr the 
model would develop roughly the observed pattern of twisting in the region between radius 
40" and 400", where our tilted-ring fit is most reliable. Because it warps away from the 
pole at larger radii, it does not become strongly twisted, as the naive arguments of 
Section [5] would suggest. 

The mass of the Hi disk is 16% of that of the stellar disk in our dynamical model. If we 
had included it in our rotation curve fit, we would have reduced the mass of both the disk 
and the dark halo to compensate. The torque from the disk would then be no more than 
16% less, and precession times would be longer by that same fraction. Our conclusion 
remains unchanged. 

Differential precession changes the twist of the gas orbits, but not their tilt with respect to 
the stellar disk. Why then should the gas disk have the observed run of tilt? In the middle 



panel of Figure [14] we show the run of tilt angle that allows the orbits at all radii to precess 
together in our dynamical model. The curve is for a halo flattened to an E3 shape, but 
those for a spherical halo and even an E6 halo lie nearby. The observed run of tilt lies 
fairly close to this curve. We conclude that the warped gas disk has the shape that it does, 
because that shape has permitted it to survive far longer than would otherwise be the case. 



The dynamical model of Figures 14 assigns the maximum plausible mass to the flattened 
disk. 



Discussion 



We have made high-resolution maps of the Hi gas in NGC 3718 and its companion 
NGC 3729. Our data cube for NGC 3718 shows multiply-peaked velocity profiles and a 
complex but highly bisymmetric structure. Using inspector, a task in gipsy, we fitted a 
tilted-ring model, in which gas following near-circular orbits about the galaxy center forms 



a warped and twisted layer. We confirm the conclusions of Schwarz (1985), that the 



prominent asymmetric dust lane marks the region where the orbits of the (dusty) gas turn 



edge-on to the line of sight. The molecular gas mapped by Krips etal. (2005) shares the 



motion of the innermost Hi gas. The unusual diffuse spiral arms fall in regions where gas 
orbits appear to crowd together on the sky. The arms are visible in blue light: new stars 



have formed in the twisted gas layer. As in other galaxies with extended Hi disks (Sancisi 



etal. 2008), spiral structure is observed far out in the disk, where self-gravity should be too 



weak to provoke the gas to instability and clumping. 



The warped and twisted Hi disk can be traced to 500" or 42 kpc from the center. It is fairly 
symmetric within 7' or 35 kpc, where the orbital period is roughly a gigayear. So the gas 
disk has probably been in place for at least a few orbits at this radius, or 2-3 Gyr. Further 
out, symmetrically-placed spiral-arm fragments to the east and west are visible in both Hi 
gas and blue light. The polar gas disk is still in the process of formation: the eastern arm 
fragment continues as a streamer of gas stretching to a cloud 60 kpc north of the galaxy 
center. Sensitive Hi maps increasingly reveal such long streamers and tails in the outer 
parts of disk systems, continuous with the regular velocity field of the galaxy, that may 



represent gas in the process of joining the galaxy (van der Hulst & Sancisi 2005). However 



the gas in polar orbit around NGC 3718 is very dusty; it is not pristine material. 

NGC 3718 has been classified as a barred Sa galaxy, but this is misleading. The apparent 
bar is an effect of looking through the edge-on disk of dusty gas, and the peculiar diffuse 
spiral arms instead represent star formation in the warped and tilted gas disk. K-band 



photometry (Tully etal. 1996) shows an exponential disk close to face-on; we find no Hi gas 



orbiting near this plane, so the old stellar disk must be almost empty of cool gas. Instead, 
NGC 3718 is typical of gas-rich early-type galaxies, where Hi gas is often found far outside 



the stellar body, and does not share the stellar kinematics (Noordermeer etal. 2005 



Morganti etal. 2006 Sparke et a/.||2008 ). When we refer our tilted-ring model for the Hi gas 
to the most probable plane for the stellar disk, the innermost gas orbits are nearly polar. 
We do not see gas orbiting in the plane of the stellar disk itself. Thus NGC 3718 is indeed 



a polar ring galaxy: as in the archetype NGC 4650A (e.g. Gallagher etal. 2002), a gas-poor 
early-type galaxy is surrounded by a highly inclined gas-rich low-surface-brightness disk. 

While the inner parts of the Hi disk are nearly polar, the outer orbits tip to lower 
inclination. This pattern of tilt minimizes the destructive effects of differential precession, 
and has allowed the polar structure to survive until the present day. The observed pattern 
of twisting can be explained by a dynamical model for the galaxy in which the gas orbits 
precess freely about the pole of the stellar disk, and the dark halo is roughly spherical. 

Polar ring galaxies are one of our few tools for studying the three-dimensional shape of the 
dark halo. Our models for NGC 3718 allows a round dark halo. St eiman-Cameron etal.\ 



(1992) obtain a similar result for the twisted dust disk in NGC 4753, concluding that 



b/a > 0.8. The Milky Way's flattening can be estimated from the near-polar streams of 
stars torn from the Sagittarius dwarf galaxy, which undergo differential precession as they 
orbit our Galaxy. This process yields confusing results, with some aspects of the streams 



pointing to a slightly oblate halo (e.g. Johnston etal. 2005) and others to a prolate halo 



(Helmi 2004): see Fellhauer etal. (2006) for a summary. However, Johnston etal. (2005) 
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favor the range 0.75 < b/a < 1.1 for the density, and strongly disfavor a very flattened halo 
with b/a < 0.6. 

By contrast, studies of the velocity fields in two polar ring galaxies imply strongly oblate 
mass distributions. When the dark halo is flattened, polar orbits are generally elongated 
towards the pole, and the polar rotation curve falls below that in the equatorial plane. By 
comparing speeds measured for gas in the polar ring with stellar speeds in the central SO 



galaxy, Sackett etal. (1994) deduced that the dark halo of NGC 4650A is considerably 



flattened, with 0.3 <b/a< 0.4, almost as flat as the stellar disk. In the system 



A0136-0801, Sackett & Pogge (1995) used Fabry-Perot imaging to map the 



two-dimensional velocity field of the polar ring in Ha emission. The kinematic major and 
minor axes were skewed away from perpendicular, a sign that the gas followed oval orbits. 
Fitting a dynamical model to the velocity field in conjunction with the spatial distribution 
of the emitting gas yielded a flattening b/a m 0.5 for the system. 

Cosmological simulations predict that the dark halos of galaxies should be fairly round. 
Halos formed from cold dark matter alone should be prolate, with axis ratios 0.6-0.7: e.g. 



Allgood etal. (2006). Adding a baryonic disk flattens the halo in the same sense as the disk 



(Dubinski 1994), but only to an average axis ratio 0.7-0.8 (Kazantzidis etal. 2004; Bailin 



etal. 2005), although Kazantzidis etal. (2004) found that 'aligned-disk' galaxy mergers 



could produce a halo as oblate as b/a ~ 0.5. If the material of the polar ring was a late 
accretion onto the central galaxy, we would expect the halo to be flattened in the same 
sense as the host's disk. If the ring gas flowed in along filaments of the 'cosmic web', as 



Maccio etal. (2006) propose, the dark halo should be oblate close to the host galaxy, 



becoming prolate and elongated along the filament further out. Bekki (1998) suggested 



that the polar ring represents the disk of a low-surface-brightness galaxy that captured the 
dense central body by merger. The dark halo might then be aligned with its long axes in 
the plane of the ring. Polar rings indeed deviate from the Tully-Fisher relation in the sense 



that Bekki's model would predict (Iodice etal. 2003); rotation speeds measured from the 



gas of the polar ring are higher than expected. 



The dark halos of galaxies may really be quite diverse; the review of Merrifield (2004) 



showed measurements covering the whole range 0.2 < b/a < 0.8. But why should we see 
such a pronounced difference among polar ring systems? Perhaps this is simply 
observational selection. Both NGC 4650A and A0136-0801 have 'classical' polar rings, 
lying nearly perpendicular to the host galaxy's stellar disk. In a galaxy with a flattened 
halo like NGC 4650A, a gas disk tipped far from the perpendicular as that in NGC 3718 
would rapidly become twisted beyond recognition. Strongly tilted rings would survive 
preferentially in systems with the roundest halos. Another possibility is that the halo 
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shape depends systematically on the galaxy's luminosity. The Milky Way, NGC 3718 and 
NGC 4753 are all luminous systems, while NGC 4650A and A0136-0801 are several times 
less luminous, with Lb ~ 4 x 1O 9 L . 
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Table 1. Basic data for NGC 3718 and NGC 3729 



Parameter NGC 3718 NGC 3729 



Type a SBa(pec) (T=l) SBa(pec) (T=2) 



Distance b 


17 Mpc 


17 Mpc 


Corrected apparent B magnitude 


10.45 


12.01 


B Luminosity 


3 x 1O 1O L 


0.7 x 10 10 L ( 


Corrected apparent K 7 magnitude 


7.35 


8.57 


K' Luminosity 


7 x 1O 1O L 


2 x 1O 1O L 


Optical radius i? 2 5 d 


3.2' 


1.3' 


Ellipticity 1 — b/a of stellar disk 


0.11 


0.32 


Inclination i { 


28° 


48° 



de Vaucouleurs etal. (1991); numerical type estimated by Tully etal. 



(1996) 



Tully (1998) 



'Extinction-corrected magnitudes M ^ and M^, from Table 2 of 



Verhei- 



jen & Sancisi (2001) 



d From D25 in Table 2 of Verheijen & Sancisi (2001) 



3 For NGC 3718, from Section 3.3 for NGC 3729, from Table 4 of Ver- 



heijen h Sancisi (2001) 



f We assume a circular disk with intrinsic axis ratio B/A = 0.2, so that 
cos 2 « = [{b/a) 2 - {B/A) 2 }/[1 - (B/A) 2 } 



Table 2. Parameters of Hi Observations 



Dates of observation 


f992.25 


R.A. pointing center (J2000.0) 


U h 32 m 39 922* 


8 pointing center (J2000.0) 


+53° 5' 25.977" 


Synthesized beam 


13.5" x 13.2" 


Effective Velocity coverage - kms -1 


708 - 1253 


Number of channels 


106 


Velocity resolution - kms -1 


5.2 / 10.4 a 


Noise - mJy/beam (rms) 


0.39 / 0.30 a 


Noise - K (rms) 


1.43 / 1.10 a 



a For the high- and low- velocity resolution cubes, re- 
spectively. 



Table 3. 20 cm Continuum Emission 



Parameter NGC 3718 NGC 3729 



Right Ascension (2000) ll h 32 m 34.934 s ll h 33 m 49.356' 
Declination 5 (2000) +53° 4' 4.63" +53° 7' 31.57" 

Flux at 20 cm (mJy) 14.4 ±1 7.9 ± 1 
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Table 4. Global Hi parameters of NGC 3718 and NGC 3729 



Parameter 


Unit 


NGC 3718 


NGC 3729 


Systemic velocity V sys a 


kms -1 


995 


1063 


Hi full width W 20 h 


kms -1 


476 


242 


Hi flux integral 


Jy kms" 1 


117.7 


3.8 


Hi radius 




8.3' = 41kpc 


1.4'= 5kpc 


Hi mass d 


1O 9 M 


8.0 


0.3 


Dynamical mass e 


1O 9 M 


500 


35 


M(Hi)/L fl 


M /L 


0.3 


0.07 



a The velocity symmetric with respect to 20% of the peak value 
on the profile wings. 

b Full width at 20% of peak flux as described in the text. 

°The total area under the global profile. 

d Calculated from the total flux assuming an optically thin 
medium. 

e For NGC 3729, adopting the inclination % from Table [TJ 
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Table 5. For the model inspector2, position angle and inclination derived for the gas 

orbits, and their tilt and twist angles relative to two reference planes: (a) that of a 
'restricted warp' and (b) that indicated by the K-band isophotes of the inner galaxy 



Radius 


PA 


i 


tilt(a) 


twist (a) 


tiltih) 


twist(b) 


arcsec 


deg 


deg 


ig = 95° 


p g = 105° 


ig = 28° 


P g = 195° 


20.0 


89.0 


93.5 


-16.1 


85.1 


100.5 


-102.6 


30.0 


98.1 


94.4 


-6.9 


84.8 


97.1 


-94.0 


40.0 


105.5 


94.8 


0.6 


111.4 


94.0 


-87.3 


50.0 


116.3 


95.0 


11.2 


89.6 


89.2 


-77.7 


60.0 


120.5 


94.9 


15.4 


89.6 


87.2 


-74.0 


70.0 


129.0 


94.6 


23.9 


90.1 


83.1 


-66.6 


80.0 


124.7 


94.8 


19.7 


89.8 


85.2 


-70.3 


90.0 


135.1 


94.1 


30.0 


90.5 


80.1 


-61.2 


100.0 


138.8 


93.7 


33.7 


90.8 


78.3 


-57.9 


120.0 


144.3 


93.1 


39.2 


91.3 


75.5 


-52.9 


140.0 


153.4 


91.7 


48.4 


92.1 


71.1 


-44.5 


160.0 


155.2 


91.4 


50.3 


92.3 


70.2 


-42.8 


180.0 


164.3 


89.7 


59.5 


93.3 


65.9 


-34.0 


200.0 


166.6 


89.2 


61.8 


93.6 


64.8 


-31.7 


220.0 


174.2 


87.5 


69.5 


94.5 


61.6 


-23.8 


240.0 


176.8 


87.0 


72.2 


94.8 


60.5 


-21.0 


260.0 


176.7 


86.7 


72.1 


95.1 


60.2 


-21.1 


280.0 


184.6 


85.0 


80.1 


96.0 


57.5 


-12.3 


300.0 


188.5 


84.0 


84.1 


96.6 


56.2 


-7.8 


320.0 


187.8 


81.1 


83.7 


99.5 


53.4 


-8.9 


340.0 


191.7 


80.1 


87.6 


100.1 


52.2 


-4.2 


360.0 


194.5 


75.2 


90.8 


104.8 


47.2 


-0.7 


380.0 


198.0 


72.8 


94.3 


106.9 


44.9 


4.0 


400.0 


193.7 


72.5 


90.3 


107.5 


44.6 


-1.8 


420.0 


193.8 


64.0 


91.1 


116.0 


36.0 


-1.8 


460.0 


197.1 


57.5 


94.4 


122.3 


29.5 


3.6 


500.0 


200.8 


49.3 


97.6 


130.5 


21.6 


12.0 
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